Quantum logic gates for coupled superconducting phase qubits 
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Based on a quantum analysis of two capacitively coupled current-biased Josephson junctions, we 
propose two fundamental two-qubit quantum logic gates. Each of these gates, when supplemented 
by single-qubit operations, is sufficient for universal quantum computation. Numerical solutions of 
the time-dependent Schrodinger equation demonstrate that these operations can be performed with 
good fidelity. 
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The current-biased Josephson junction is an easily fab- 
ricated device with great promise as a scalable solid-state 
qubit Q, as demonstrated by the recent observations of 
Rabi oscillations This phase qubit is controlled 

through manipulation of the bias currents and applica- 
tion of microwave pulses resonant with the energy level 
splitting 0- 

In this Letter we analyze the quantum dynamics of two 
coupled phase qubits. (The classical dynamics of this sys- 
tem has also been studied recently |^]). We identify two 
quantum logic gates that, together with single-qubit op- 
erations, provide all necessary ingredients for a universal 
quantum computer. We perform full dynamical simula- 
tions of these gates through numerical integration of the 
time-dependent Schrodinger equation. These two-qubit 
operations may be experimentally probed with the meth- 
ods already used to observe single junction Rabi oscilla- 
tions |3, H|. Such experiments are of fundamental im- 
portance: the successful demonstration of macroscopic 
quantum entanglement holds profound implications for 
the universal validity of quantum mechanics 0- Im- 
portant progress toward this goal are the temporal os- 
cillations of coupled charge qubits and spectroscopic 
measurements |]| on the system considered here. Finally, 
our methods are applicable to the other promising super- 
conducting proposals based on charge, flux, and hybrid 
realizations |8|. 

Figure ^a) shows the circuit diagram of our coupled 
qubits. Each junction has characteristic capacitance Cj 
and critical current I c , and they are coupled by capaci- 
tance Cc ■ The two degrees of freedom of this system are 
the phase differences _Tl and 72 , with dynamics governed 
by the Hamiltonian [9( 



H 



4E C (1 + Q- l h- 2 (p{ +p§ + 2Cpip 2 ) 
-i?,/(cos7i + J171 + cos 72 + J27 2 )- 



(1) 



Here we have employed the charging and Josephson ener- 
gies Ec = e 2 /2Cj and Ej = hl c /2e, the normalized bias 
currents J\ — I\/I c , Ji = Iijlc, and the dimensionlcss 
coupling parameter £ = Cc/{Cc + Cj). 

This coupling scheme has been recently analyzed [f^liol 
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FIG. 1: Capacitively coupled Josephson junctions: (a) ideal 
circuit diagram, (possible experimental parameters include 
Cj = 6 pF, C c = 60.6 fF, and I c = 21/iA, J = 0.988, 
u}o/2n = 6.45 GHz); (b) time-dependent ramp of bias cur- 
rents, specified through the detuning e (see text). 



and results in a system with easily tuned energy lev- 
els and adjustable effective coupling. While £ is typically 
fixed by fabrication, the energy levels and the effective 
coupling of the associated eigenstates are under experi- 
mental control through J\ and J2- As shown below, the 
two junctions are decoupled for J± and J \ sufficiently 
different, but if J\ and Ji are related in certain ways, 
the junctions are maximally coupled. To illustrate this 
method of control, we define a reference bias current Jo 
and consider the variation of J\ and J 2 through a detun- 
ing parameter e: 



VT - J2= Vl — Jo(l-e). 



(2) 



Quantum logic gates are implemented by varying e 
with time as shown in Fig. ^b). This ramps the bias 
currents, moving the system smoothly (with ramp time 
tr) from 6a, where the eigenstates are essentially un- 
cntangled, to eg, where the eigenstates are maximally 
entangled. Entangling evolution is then allowed to oc- 
cur for an interaction time 17, after which the system is 
ramped back to 6a- 

For analysis, we use the energy scale hujQ = 
■\/8EcEj(l — Jo) 1 / 4 ( w o/27r is the classical plasma fre- 
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quency of a single junction) and the effective number N s 
of single junction (metastable) energy levels 



(3) 



After choosing a fixed coupling £, gate design requires 
the identification of suitable N s , e A , €b, Tr, and 17. Fig- 
ure |2] shows the relevant energy levels of H as a func- 
tion of e for the physically interesting case C = 0.01 and 
N s = 4, and with the potential energy minimum sub- 
tracted off. The energy levels E n (e) and their associated 
two-junction eigenstates |n;e) were computed using the 
method of complex scaling || 0] applied to the cubic 
approximation [13l | of H. 

In general, each energy state \n; e) is an entangled su- 
perposition of the product states \jk; e) = |j; e) ® |/c; — e), 
where |j; e) are energy states of an isolated junction with 
normalized bias current J\ . (The "round" and "angular" 
brackets distinguish the coupled and uncoupled bases, 
respectively.) However, for |e| > 0.1, the energy states 
are essentially unentangled and well approximated by the 
product states, which are used to label the corresponding 
energy levels in Fig. [21 Thus for e A = —0.1 we find that 
the eigenstates satisfy the relations |1;£.a) — |10;£a)) 
\2;e A ) = \01;e A ), and \A;e A ) S |ll;e A ). The ground 
state |0;e) = |00;e), not shown, is essentially unen- 
tangled for all e. We choose these states for our two- 
qubit basis. In addition, there are the auxiliary states 
|3;e A ) £* \20;e A ) and \5;e A ) £* |02;e A ). 

For e near e± = ±0.04 and 0, where avoided level 
crossings occur, we find significant entanglement. Fig- 
ure [S] shows the entanglement of the states \n;e), with 
n = 1,3,4,5, as a function of e. (The entanglement of 
states 1 and 2 are nearly identical.) The entanglement is 
given in ebits ^ij: a state with one ebit entanglement is 
a maximally entangled two-qubit state. The gates con- 
structed below use this entanglement to perform two- 
qubit operations. 

It is important to account for the nonqubit states |02; e) 
and 1 20; e) — poorly designed interactions will result in un- 
wanted evolution of 1 11; e) into these auxiliary levels, an 
effect called leakage. However, we show below that at 
both cb = and es = £- the state |ll;es) is a superpo- 
sition of only two of the energy eigenstates. Therefore, its 
time evolution is of the form |(11; eB\e~ lHT '^ h \ll; eg) | 2 = 
a + b cos 2 (Qti), with a + b = 1. Choosing tj = kir/Q,, 
where k is an integer, ensures that the oscillation of 1 11; e) 
with the auxiliary states completes k full cycles, mini- 
mizing leakage. This procedure is similar to operations 
performed in ion traps fl5| . 

If we allow 1 1 1 ; e) to evolve through states |02;e) and 
1 20; e), we must consider another possibility for error: the 
tunneling rates of these auxiliary states are higher than 
that of 1 11; e). As these tunneling rates are all proportial 
to e ~ 36N =/ 5 W e minimize this error by choosing a 
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FIG. 2: Normalized energy levels, N s — 4, £ = 0.01, as a 
function of detuning parameter e: (a) the energy levels E\ (e) 
and .E2 (e) with avoided level crossing at e = 0; (b) the energy 
levels i?3(e) through £7s(e) with avoided crossings of 4 and 5 
at e± ~ ±0.04, and 3 and 4 at e = 0. The ground state en- 
ergy Eo (not shown) is approximately 0.981/icjo- The natural 
two-qubit states of this system are |0; e A ), |1; ea), |2; e A ) and 
|4; e A ), with e A = — 0.1 (see text). 
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FIG. 3: Entanglement, N 3 = 4, f = 0.01: The entanglement 
of the energy eigenstates as a function of detuning parameter 
e: (a) |l;e); (b) |3;e); (c) |4;e); and (d) |5;e). 



large N s . For example, to achieve a fidelity greater than 
0.999, we find that N s > 4 is necessary. We note that 
proper operation requires the ramp time tr to be adi- 
abatic with respect to single junction energy level spac- 
ings, yet nonadiabatic with respect to the coupled energy 
level splittings. This leads to our choice C = 0.01. 

Before analyzing gate dynamics, we look more closely 
at the eigenstates at e = and e±. These can be un- 
derstood perturbatively and identified with degeneracies 
in the combined spectrum of two uncoupled Josephson 
junctions. For example, with £ = at e = the states 
|01;0) and [10; 0) are degenerate. Coupling splits this 
degeneracy [lfj, a,nd yields the true eigenstates 



|1;0) = 2- 1 /3(|01;0)-|10;0» 
|2;0) = 2^/2(101; 0) + |10;0)). 



(4) 
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These states (Fig. [31 a)) are maximally entangled qubits. 
Similarly, states |02;0) and |20;0) are degenerate with 
C = 0. With nonzero coupling, however, the energy split- 
ting in this case is much smaller [see the inset to Fig|5Jb)]. 
Here, a three-state analysis is required, with the nearly 
degenerate |11;0) mediating the coupling. The eigen- 
states are approximately 

|3;0) = 2- 1 / 2 cos6>(|02;0) + |20; 0» - sin0|ll; 0) 

|4;0) = 2-V2(|02;0>-|20;0» (5) 
|5;0) = 2- 1 / 2 sin6l(|02;0) + |20; 0)) + cos0|ll; 0), 

with 9 = 0.185. In accord with Fig.EJb-dh each of these 
states is substantially entangled at e = while |11; 0) 
is a superposition of |3;0) and |5;0). 

Finally, perturbation theory shows that the off- 
symmetry degeneracies occur for e± = ±5/367V s . For 
example, at e = e_ with £ = states |02; e_) and | XI; e_) 
are degenerate. Coupling leads to the true eigenstates 



E 2 (0) - 2E (0)]t i /2h. For U 2 we find the swaplike gate 



|4; e _) = 2-V2(|02;e_)-|ll; e _)) 
|5; e _) = 2- 1 /2(|02;e_) + |ll;e_)). 



(6) 



As seen in Fig. EJc,d), both eigenstates have an entan- 
glement of 1 ebit. Further, we see that |ll;e_) is a su- 
perposition of |4;e_) and |5;e_). 

We now describe the two-qubit gate that uses the en- 
tanglement at e_. As noted above, if the state 1 11; e_) 
is prepared, its time evolution (for fixed e) will be 
oscillatory. Letting tj = 2irh/[Ec ) (e^) — £^(e_)] = 
V^vr/Cwo, |ll;e_) performs a complete oscillation, while 
picking up an overall controlled phase. The remain- 
ing states also evolve dynamical phases, which can 
be factored out as one-qubit gates by letting U\ — 
e lai R z {a 2 )® R z (a 3 )e- lHr '/ h , with R z (9) = e - i6a */ 2 (a z 
is a Pauli matrix). Here a\ — \E\(eJ) + E 2 (c-)]ti /2H, 
a 2 = [Ei(e-) - Eo(e-)]Ti/h, and a 3 = [E 2 (e-) - 
Eq(c— )]ti/K. In our two-qubit basis, this operation is 
the controlled-phase gate 



/l 00 

10 

1 

\ e~ z4 



J 



(7) 



with (j) = [E 4 (e-) + Bo(e-) - E x (eJ) - E 2 {e-)]T I /h. For 
N s = 4 and C = 0.01, we find <fi — 1.027T, thus this gate 
is approximately the controlled- Z gate |l6| . 

We can also use the entanglement at = for 
quantum logic. From the dynamics of 1 11; 0) [im- 
plied by Eq. J5J], we let the interaction time be tj = 
2irkh/[Er,(0) — E 3 (0)}, where k is an integer. From 
Eq. (0J, however, the states |01;0) and 1 10; 0) will also 
oscillate 0. Removing one-qubit dynamical phases as 
above, we define U 2 = e lai R z (a 2 ) ® R z (a 3 ) e - lHTl l h , 
with ai = [^i(0) + E 2 (0)}t i /2H and a 2 = a 3 = [E^O) + 
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(8) 



with swap angle 9\ = [£^(0) — Ei(0)]ti/2H and 
controlled-phase 9 2 = [^ 5 (0)+^o(0)-^i(0)-B 2 (0)]Tj/n. 
As 9\ and 9 2 are in general irrational multiples of it, 
this gate is universal for quantum computation |l7j . For 
example, by tuning Jo such that N s = 5.16 and let- 
ting k = 2, the full swap dynamics is generated, with 
6*i = 7r/2, 9 2 = 7r/4, and tj = k/Quj q . 

Proceeding beyond this heuristic analysis, we have nu- 
merically solved the time-dependent Schrodinger equa- 
tion using split-operator unitary integration 18]. Tun- 
neling is incorporated through an absorbing boundary 
condition 19]. Taking = —0.1 as our initial detun- 
ing, we evolved states having initial conditions |0;ca), 
|1; €a), |2; e^), and |4; ca) using the the ramp function of 
Fig.[TJb). While the results quoted below are in the cubic 
approximation, results obtained using the full Hamilto- 
nian are only marginally different. 

The controlled-phase gate U\ is simulated with N s — 4 
and ( — 0.01. Using the minimum splitting between 
Ei(<e) and £s(e) we take es = —0.036, with tr — 20tt/loo 
and ti — 434/o7o- The dynamical behavior of this gate 
is illustrated in Fig. Ufa), where the probability P(t) = 
\{A]e A \^{t))\ 2 with |V>(0)) = \^ £ a) is shown. Including 
all two-qubit states, we find an average gate fidelity 20] 
F = 0.996, with a leakage [5l[ L = 0.003. 

For the swaplike gate U 2 , we use _/V s = 5.16, £ = 0.01, 
tr = 207r/o;o and tj — 278/wo. The probability P(t) is 
shown in Fig. Etb), completing two oscillations (k = 2) 
but with diminished amplitude due to the shift of the po- 
tential (since |(11; e^|ll; e B )\ 2 < 1). The swap of 1 10; €a) 
and [01; ca) is shown in Fig. This gate is not quite 
as good as U\, with a fidelity F = 0.972 and leakage 
L = 0.006. We expect that a complete optimization of 
U\ and U 2 will yield sufficient fidelity for fault-tolerant 
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FIG. 4: Dynamical evolution of state with initial condition 
|4; ea). The probability P(t) = |(4; e A \ip(t))\ 2 is shown for (a) 
the phase gate Ui and (b) the swaplike gate U2- 
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FIG. 5: Dynamical evolution of state with initial condi- 
tion |1;6a) under the swap like gate. The contours represent 
the numerically computed wave function (modulus squared) 
evolving from what is nearly 1 10; ca) to |01; 6a). 



qubit basis generates a controlled phase on state 1 1 1 ; e) , 
while for U2 an additional swap operation is performed 
between states 1 01 ; e) and 1 10; e). Finally, we expect that 
a 3-junction design, using the unitary operations identi- 
fied here, is an even better candidate for universal quan- 
tum computation with capacitively coupled Josephson 
junctions. 

We would like to thank A.J. Berkley, J.M. Martinis, 
R. Ramos, H. Xu, and M. Gubrud for useful discussions. 
This work was supported in part by the U.S. Department 
of Defense and the State of Maryland through the Center 
for Superconductivity Research. 



quantum computation 22] . 

An experimental demonstration of these gates requires 
the following considerations. First, for typical experi- 
mental parameters [see Fig. ^a)] the plasma frequency 
ujq/2tt is near 6 GHz. With this value, the ramp time 
used above is tr = 1.67 ns. The total gate times 
(tj + 2tjj) are t\ = 14.85 ns for the phase gate and 
t 2 = 10.7 ns for the swap gate. Pulse shapes similar to 
Fig. ^b) with these times can be engineered with con- 
vential electronics. Second, coherent operation at these 
time scales requires low dissipation due to the control 
circuit, which is possible with impedance transformers 
• Finally, the essential parameters controlling the gate 
dynamics are the energy level spacings. Since these can 
be determined spectroscopically 0, Ej, all aspects of this 
design are experimentally accessible. 

Other important issues are noise in the bias currents 
and nonidentical junction parameters. Current noise will 
cause fluctuations in e, while nonidentical junctions will 
have reduced symmetry about e = 0. The TJ% gate is 
particularly sensitive to these, as it uses both the sym- 
metry and delicate structure of the eigenstates at e = 0. 
In contrast, since the phase gate U\ operates at an off- 
symmetry position, its operation will be less sensitive to 
these sources of decoherence. 

The controlled coupling of qubits can be refined by 
introducing a middle junction to generate entanglement 
between adjacent qubits [ljj. Then, on a state of the 
form I*) = c 00 |000) + c i|001) + ci |100) + c n |101), the 
operation U = {U 2 ® I) (I ® Ui)(U 2 ® /) (where / is the 
one-qubit identity operator) is equivalent (with single- 
qubit operations on the outer junctions) to performing a 
phase gate on the outermost qubits, leaving the central 
qubit in its ground state. In this scheme the swaplike 
gate U2 does not act on 1 1 1 ; e) , so decoherence effects 
will be less severe than indicated above. Furthermore, 
a larger coupling may be used, leading to smaller gate 
times. 

In conclusion, we have shown how to implement two 
quantum logic gates in this coupled junction system. For 
Ui , evolution through auxiliary levels outside of the two- 
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